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ABSTRACT 



Context. The properties of the X-ray emission lines are a fundamental tool for studying the nature of the matter surrounding the 
neutron star and the phenomena that produce these lines. 

Aims. The aim of this work is to analyze the X-ray spectrum of 4U 1538-52 obtained by the XMM-Newton observatory and to look 
for the presence of diagnostic lines in the energy range 0.3-1 1.5 keV. 

Methods. We used a 54 ks PN &■ MOS/XMM-Newton observation of the high mass X-ray binary 4U 1538-52 covering the orbital 
phase between 0.75 to 1.00 (the eclipse-ingress). We have modelled the 0.3-11.5 keV continuum emission with three absorbed power 
laws and looked for the emission lines. 

Results. We found previously unreported recombination lines, in this system, at ~2.4 keV, ~1.9 keV and ~1.3 keV, consistent with 
the presence of highly ionized states of S XV Hea, Si XIII Heo- and Mg Ka or Mg XI Hea. On the other hand, both out of eclipse 
and in eclipse we detect a fluorescence iron emission line at 6.4 keV which is resolved into two components: a narrow (cr < 10 eV) 
fluorescence Fe Ka line plus one hot line from highly photoionized Fe XXV. 

Conclusions. The detection of new recombination lines during eclipse-ingress in 4U 1538-52 indicates that there is an extended 
ionized region surrounding the neutron star. 
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1. Introduction 

4U 1538-52, discovered by the Uhuru satellite 
jGiacconi et a l. 1974), is an X-ray pulsar with a B-type 
supergiant companion, QV Nor. It has an orbital period of 
-3.728 days dPavison et al. 19771 ICTark 2000l l. with eclipses 
lasting ~0.6 days dBecker et al. 197'7] l. This X-ray persistent 
system produces this radiation when the neutron star captures 
matter from the wind of the B supergiant star. Assuming a dis- 
tance to the source of 5.5 kpc (Becker etal. 1977J and isotropic 
emission, the estimated X-ray luminosity is ~ (2 - 7) X 10^^ erg 



In 4U 1538-52 the emission line at ~6.4 keV can usually 
be described within the uncertainties either by a single nar- 
row Gaussian line or by a multiplet of narrow Gaussian lines 
(W hite et al. 1 983). Observations carried out with Tenma de- 
tected it at 6.3+0.2 keV and EW 50+30 eV (Makishima et all 
1987), while the Rossi X-ray Timing Explorer (RXTE) saw it at 
6.25+0.06 keV and EW 61 eV (Coburn et al. 2002). Other X-ray 
observatories such as BeppoSAX (Robba et al. 2001 ) and RXTE 
have used a Gaussian line at ~6.4 keV for describing the flu- 



orescence of iron in a low-ionization region (Mukherjee et al. 
2006; IRodes 20071 ). The variability of this line was studied by 
^ ■ s~' in the 3-100 keV range (IRodes 20071 1. Therefore, the size |Rodes etal. (2006)| 



of the ionization zone may be a relatively small region in the 
stellar wind ( [Hatchett & McCray 1977[lvan Loon et al. 20011) . 

Fluorescence iron emission lines from X-ray pulsars are pro- 
duced by illumination of neutral or partially ionized material by 
X-ray photons with energies above the line excitation energy. 
Possible sites of fluorescence emission may be: (i) accretion disk 
(mostly seen in low mass X-ray binary sources); (ii) stellar wind 
(in high mass X-ray binary pulsars); (iii) material in the form 
of a circumstellar shell; (iv) accretion column; (v) material in 
the line of sight; or some combination of these locations. In this 
sense, fluorescence lines in the X-ray spectrum are an interest- 
ing tool for studying the surrounding wind regions and elemental 
abundance in X-ray sources. 
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In this paper, we present a spectral analysis based on the ob- 
servation of 4U 1538-52 performed with data from the XMM- 
Newton satellite. The observation covers the orbital phase inter- 
val 0.75-1.00 and we detect the presence of the K^, iron line 
at ~6.4 keV and some blended emission lines below 3 keV. In 
Sect.|2]we describe the observation and data analysis. We present 
in Sect. [3] timing analysis; in Sect. |4] spectral analysis, and in 
Sect. |5] we summarize our results. 

2. Observation 

4U 1538-52 was observed with the XMM-Newton sateUite in 
2003 from August 14 15:34:01 to August 15 14:02:30 UT using 
EPIC/PN and EPIC/MOS. The three European Photon Imaging 
Cameras (EPIC) consist of one PN-type CCD camera and 
two MOS CCD cameras dStruder et al. 200UlTumer et al. 20011 . 
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Table 1. Log of XMM-Newton observation of 4U 1538-52. 



Instrument 


MJD 


Exposure 


Orbital 




Start-End 


(ks) 


phase 


PN 


52 865.66-52 866.58 


53.68 


0.75-1.00 


MOS-1 


52 865.65-52 866.58 


64.33 


0.75-1.00 


MOS-2 


52 865.65-52 866.58 


67.93 


0.75-1.00 



3. Timing analysis 

In Fig. [1] we show the background subtracted light curves for 
the EPIC/PN camera, binned at 20 s, in the energy ranges 0.3-3 
keV and 3-11.5 keV and their hardness-ratio (HR) calculated 
as (3-1 1.5 keV)/(0.3-3 keV) (hard/soft). This light curve shows 
the source out of eclipse during the first ~25 ks and in eclipse 
in the last ~55 ks. The shape of the light curve is similar both in 
the hard and soft components during the eclipse. During the first 
~10 ks there is a softening trend but the HR presents significant 
differences between ~10 ks and ~20 ks. In this time interval, 
the soft component has a count rate similar to that in eclipse, 
whereas the hard component appears uneclipsed. 

The XMM-Newton observation of 4U 1538-52 allows us to 
investigate the pulse profile in detail. To estimate the pulse pe- 
riod, we used the XRONOS version 5.21 timing analysis soft- 
ware package. First of all, we estimated the pulse period with 



the powspec task; secondly, we searched for the best pulse period 
with the efsearch task using the test; and finally, we obtained 
the folded light curve at the best-fit period P - 526.7 + 0.2 s. 
This result is consistent with the RXTE pulse period obtained by 
Mukherjee et"ai:|(2006), P = 526.849 + 0.003 s. Fig.|2]shows the 



The three EPIC cameras were operated in Full frame mode, with 
a time resolution of 73.4 ms for the PN camera and of 2.6 s for 
the two MOS cameras, and the thin filter 1 {MOS-1 and PN) and 
medium filter {MOS-2) were used. The net count rate of each one > 
was 3.1 16 ±0.008 ctss-i, 0.723 + 0.003 ctss-i and 0.734 ±0.003 7 
cts s"', respectively. The reflection grating spectrometer {RGS) 3 
showed a low level of counts, 1.4 ± 1.1 x lO^^* cts s"', and we 
did not use this spectrum in our analysis. The resulting effective 
exposure time was 53.68 ks for the PN, 64.33 ks for the MOS-1 | 
and 67.93 ks for the MOS-2. The details of the observation are 12 
listed in Tabl^U To calculate the orbital phase we used the best- J 
fit ephemeris of lMakishima et aL] ("1987) with the orbital period 
Porh = 3.72854 ±0.00015 days and the eclipse centre 45,518.14 
days (Modified Julian Date). - 

We reduced the data using Science Analysis System (SAS) J 
version 8.0, using the most up-to-date calibration files. Spectra ^ 
and light curves were extracted using circular regions centered 
on the source position, with radii of 600 pixels (30 arcsec) for 
EPIC/PN and 400 pixels (20 arcsec) for the two EPIC/MOS in- 
struments. We selected background data from circular regions 
off'set from the source, with radii six times the source extraction 
radius. We accumulated all the events with patterns 0^ and 0- 
12 for EPIC/PN and EPIC/MOS, respectively. 

As 4U 1538-52 is a bright source, we analyzed whether our 
data were affected by pile-up, i.e. if more than one X-ray photon 
arrives at either a single CCD pixel, or adjacent pixels, before the 
charge is read out. This effect can distort the source spectrum and 
the measured flux. The XMMSAS provides the task EPATPLOT 
to determine the distribution patterns produced on the CCDs by 
the incoming photons and verify whether the source is piled-up 
or not. We created an event list extracted for both an annulus and 
a circle region and displayed the fraction of single, double, triple 
and quadruple pixel events. The resulting plots of the PN and 
MOS data revealed no obvious pile-up, even when using a circle. 
The maximum count rate before pile-up is typically ~0.7 cps 
for the MOS instruments and 6 cps for PN in full frame mode. 
Therefore, we conclude that our data are not affected by pile-up. 



folded light curves in the two energy intervals 0.3-3 keV (soft) 
and 3-11.5 keV (hard) together with the folded hardness-ratio. 
The HR shows two peaks at pulse phases 0.1 and 0.6, respec- 
tively. 




Fig. 2. EPIC/PN folded fight curves of 4U 1538-52, binned at 
20 s, in the energy ranges 0.3-3 keV (top panel), 3-11.5 keV 
(middle panel) and their hardness ratio (3-1 1 .5)/(0.3-3) (bottom 
panel). 



4. Spectral analysis 

For spectral analysis we used the XSPEC version 12.5.0 
(Arnaud 1996) fitting package, released as a part of XANADU 
in the HEASoft tools. 

The EPIC cameras cover the energy range 0.3-11.5 keV. 
In previous works based on a wider energy range {Tenma, 
1.5- 35 keV, IMakishimaetaT] 1987; Beppo SAX, 0.1-100 
keyiRobbaetal. 2 001; RXTE, 3-100 keV, Coburn et al. 20021 
IRodes 20071 INTEGRAL, 3-100 keV, iRodes-R oca et al. 200§ 
the X-ray continuum of 4U 1538-52 was described by differ- 
ent absorbed power law relations, modified by a high-energy 
cutoff plus a Gaussian emission line at ~6.4 keV, when- 
ever present. Moreover, these models have also been mod- 
ified by absorption cyclotron resonant scattering features at 
-21 keV dCTark et al. 19901 IRobba et al.^OOTT) and at -47 keV 
jRodes-Roca et al. 20091 ). 

Other X-ray continuum models of high-mass X-ray binaries 
(HMXBs), such as two absorbed power laws or three absorbed 
power laws, can also describe the observations with ASCA of 
Vela X-1 jSako et al. l"999]l or wifli Chandra & XMM-N ewton 
of 4U 1700-37 (Boros on et alTIOOSHvander Meer et al.l 2005). 
One power law is used to fit the direct continuum which origi- 
nates near the neutron star Part of this radiation is scattered by 
the stellar wind of the supergiant star and is fitted by the sec- 
ond power law. Some HMXBs present a soft excess at low en- 
ergies (0.1-1 keV) which may be modelled as another absorbed 
power law, blackbody component orbremsstrahlung component. 
The physical origin of this component is not well understood, al- 
though it is reported for some X-ray binary sources, such as 4U 
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Fig. 1. EPIC/PN light curve of 4U 1538-52, binned at 20 s, in the energy ranges 0.3-3 keV (top panel) and 3-1 1.5 keV (middle 
panel), and their hardness ratio (bottom panel). 



1700-37 jvan der Meer"etari 2005), 4U 1626-67, Cen X-3 or 
Vela X-1 dHickox et al. 20041) . The BeppoSAX spectrum of 4U 
1538-52 showed a soft excess at 0.3-1.0 keV, which was mod- 
elled with a blackbody component with a temperature of 0.08 
keV dRobba e t al. 2001 ). Moreover, they also obtained a good fit 
using a bremsstrahlung component with a temperature of 0.10 
keV. We used all of these models to fit the XMM-Newton con- 
tinuum spectrum of 4U 1538-52 and Gaussian functions to de- 
scribe the emission lines. 

4. 1 . X-ray continuum 

In order to look for diagnostic emission lines a proper de- 
scription of the underlying continuum is critical. Because of 
its higher count rate we used only the PN spectrum to pin 
down the parameters of the phase averaged spectrum and 
then we combined both PN and MOS pulse phase averaged 
spectra in order to check our results. We fitted the 0.3-1 1.5 
keV energy spectrum using three absorbed power-law compo- 
nents fixing the photon index parameter to the same value for 
each power law component, but with different normalization 
and column density Nh- We modelled the hard X-ray continuum 
(3-11.5 keV) using two power laws, each of them modified by 
different absorbing columns dMorrison & McCammon 19831 1. 



We obtained adequate fits using the same photon index a ~ 
1.134^Q Q23 for both components which describe the X-ray con- 
tinuum radiation from the neutron star: a direct emission from 
the compact object, absorbed through the surrounding stellar 
wind, and scattered radiation produced by Thomson scattering 
by electrons in the extended stellar wind of the supergiant star 
Taking the hydrogen columns listed in Table |2] into account, 
the hard component is assumed to originate from the near 
surroundings of the neutron star and is highly absorbed by 
optically thick matter in the line of sight {Nh ~ 6.0 x 10^^ 
cm^), and the scattered component is produced by optically 
thin plasma (A^H ~ 1.2xl02^cm 2).OtherHMXBswithspec- 
tra_which^anbe described by these components are Cen X- 
3 (INagase et al.| 1992; Ebisawa et al. 19961 and 4U 1700-37 
(van der Meer et al. 2005). These systems also showed a hard 
component highly absorbed even during the eclipse. 

We used the third absorbed power-law because a soft excess 
at 0.3-1 .0 keV is detected in the PN spectrum. We fixed the pho- 
ton index to that of the other power-laws obtaining a significant 
improvement (xl - 1.16, with an F-test of 4.2x10"'""). Fig. |3] 
shows the three different absorbed power-law continuum com- 
ponents and their residuals in units of cr. The spectrum clearly 
shows several emission lines at 6.4 keV, 2.4 keV, 1.9 keV and 
1.3 keV, and an absorption feature at 2.1 keV. The flux ratios 
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Table 2. Fitted parameters for the PN X-ray continuum spectrum 
in Fig. [3] 



Component 


Parameter 


3PL'' 


2PL+bb* 


Hard 


Nh (10^- cm-^) 
Photon index^ 
Normalization^ 


60 ±4 

1 1^4+0.018 
^•^•^ -0.023 

93 ± 9 


62+* 

1 1 O+0.03 
^•^-'-0,05 

92+1' 


Scattered 


Nh (10^' cm"^) 
Normalization' 


ii.70+;;:j' 

43 ± 3 


12.7!»;^ 
47+f 


Soft 


Nh (10^2 cm-2) 
Normalization'' 
kTtb (keV) 
Normalization;,;, 


0.74 ± 0.04 
7 4.7+0.15 


42+"-'''' 

925 +"'"5 
1.8+ij X 10-^ 




^J(dof) 


1.16(2235) 


1.16(2234) 



" Three absorbed power-laws 

* Two absorbed power-laws plus an absorbed blackbody 
The same for all the power-laws 
(x 10"*) in units of photons s"' cm"' keV"' at 1 keV 



are usually above 30%. We are sure that the low energy residu- 
als are not due to incomplete calibration. The cross-calibration 
XMM-Newton database consists of ~150 observation of differ- 
ent sources, optimally reduced, fit with spectral models defined 
on a source-by-source basifl As a result, deviations for PN and 
MOS flux ratios above 10% are consistent with the presence of 
emission or absorption lines. 



0.1 





/ 

/ 


: 1 "1 1 1 1 1 ^ 









0.5 I 2 5 10 

channel energy (keV) 

Fig. 3. EPIC/PN spectrum of 4U 1538-52 and X-ray continuum 
model. The dash line represents the direct component, the dash- 
dot-dot-dot the scattered component, and the dot line the soft 
component (top panel). The lower panel shows the residuals be- 
tween the spectrum and the model. 



We also used other components to investigate the physical 
origin of this soft-excess, such as a bremsstrahlung or blackbody 
component. These model components need to be modified by an 
absorption column too. Although a bremsstrahlung component 
improves the fit significantly (xl - 1.17, F-test of 4.0x10"'""), 
its inferred temperature is too high (196.9 keV) and could not be 



constrained within the XMM band pass, suggesting that the soft 
excess has another physical origin. 

In Table |2] we report the model parameters of the X-ray con- 
tinuum spectrum for the PN camera. All uncertainties refer to a 
single parameter at the 90% (Ax^ - 2.71) confidence limifl. 

When we used a blackbody component we also obtained a 
similar improvement of the fit (xl = 1.16, F-test of 1.6x10"'^^). 
The blackbody has a temperature of kT ~ 0.925+° °[^ keV. We 
calculated the luminosity of the source of the blackbody compo- 
nent over the 0.3-1 1.5 keV energy band, assuming a distance to 
the source of 5.5 kpc jBecker et al. 197'7b and isotropic emission 
to be 5.4X10''-' erg s"'. Since the surface luminosity of a black- 
body depends only on its temperature, it is possible to calculate 
the size of the emitting region: 



Rhbikm) = 3.04 X lO'* 



^ bh 



(1) 



where D is the distance to the source in kpc, Fht the unabsorbed 
flux over the 0.3-1 1.5 keV energy band in erg s"' cm"^, T^i, the 
temperature in keV and Rhb is the radius of the emitting region. 
From Eqll] we found a radius of the emitting surface of 0.24 km. 
If we assume thermal emission from the neutron star polar cap, 
this radius may be consistent with the expected size. However, 
we estimated the radius of the accreting polar cap using the equa- 
tion (Hickox et al. 2004) : 



RrnI — Ra 



Rm 



(2) 



where /?,„ is the magnetospheric radius and R^s the neu- 
tron star radius. The magnetospheric radius may be estimated 
from the magnetic dipole moment and mass accretion rate 
(lElsner & Lamb 19771 [Audley et al. 1996) 1: 



R,n(cm) 



3.2 X 10** M^7^^ ull^ 



M 

Mo 



-1/7 



(3) 



where fi^o is the magnetic dipole moment of the neutron star in 
units of 10^" G cm^ and Mu is the mass accretion rate in units 
of 10'^ g s"'. With M ^ 1 3Me (Reynol ds et al. 1992|l, M n = 
883 g s"i dClark et al. l"99T;'Rode s 2007iiRodes et al. 20081 1 and 
yL<3o = 1.8 G cm^ (Rodes 20071), the Eq. |3] yields the value 
R,,, ~ 6.2x10^ cm. Assuming R^s ~ 10^ cm, we obtained 
Rcoi ~ 1 .27 km, five times our estimated blackbody emitting ra- 
dius. Although this result suggests that the soft excess is not 
formed by blackbody emission, taking into account the asso- 
ciated errors in the parameters, we can not exclude black- 
body emission as the origin of the soft excess. 

We also tried to fit this soft excess by adding 2 to the photon 
index of the soft component compared to the index of the hard 
component, i.e. scattering by dust grains at large distance from 
the source (IRobba et al. 2001llvan der Meer etal., 2005). We ob- 
tained a good description of the spectrum (xl = 1 -20, F-test of 
11x10""^), but a significantly worse fit than the previous com- 
ponents (xl = 1.16). 

Finally, we modelled the soft X-ray continuum compo- 
nent with a blend of Gaussians only (IBoroson et al. 20031 
Ivan der Meer et alTl 2005) and refitted the data. We could not ob- 
tain a good description below 0.6 keV, though. Therefore, al- 
though the model composed of three power-law components 
also does not fit the data well below 0.6 keV, it was the best 
model obtained and we used it for our analysis. 



' http://www.iachec.org/meetings/2009/Guainazzi_2.pdf 
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Afterwards we combined EPIC/PN and EPIC/MOS spec- 
tra and applied the same models to them obtaining similar 
results. The spectra of all three EPIC instruments were fit- 
ted simultaneously, including a factor to allow for the adjust- 
ment of efficiencies between different instruments. We fixed 
the continuum parameters, described in Tabled and looked 
for the emission features (see Sect. 14.2) . 

4.2. Emission lines 

The X-ray eclipse spectra of some HMXBs show many emis- 
sion lines and radiative recombination continua (e.g. Vela X- 
IJSakoetal. 1999. Cen X-3 Ebisawa et al. 199 6: 4U 1700-37 
Ivan der Meer et alT] 2005). The data together with the best-fit 
model, three absorbed power laws plus emission lines, and resid- 
uals of the fit as the difference between observed flux and model 
flux divided by the uncertainty of the observed flux, i.e. in units 
of cr, are included in Fig. H] As we can see in Fig. H) the PN 
spectrum of 4U 1538-52 shows emission lines below 3 keV 
and between 6 and 7 keV. We fitted these lines as Gaussian 
profiles (see Table [3] for parameter values). Because the energy 
resolution of the EPIC cameras is not sufficient to resolve all 
lines (FWHM 80-100 eV in the energy range 1-3 keV), many 
of them could be blended (e.g., 1.75 keV Si Ka, Al xiii Lya 
and 1.85 keV Si xiii Hea). Using the list of emission lines in 
Ivan der Meer et alTl (2005). we could identify suitable candidates 
among fluorescence emission lines from near-neutral species 
and discrete recombination lines from He- and H-like species 
(see also Dra ke 1988l and |Hojnacki et al. 2007| . 

Fig. ID clearly shows the presence of the Fe Ka line at ~6.4 
keV, but a second broad emission line around 6.65 keV or an ab- 
sorption edge at ~7.2 keV associated with low-ionization iron is 
an unresolved issue within the deviation for PN and MOS ratios. 

We started by modeling the fluorescence iron line (see Fig.|4]l 
that was also detected by several other satellite observatories, 
such as Tenma, BeppoSAX or RXTE. Starting with the 3PL con- 
tinuum model, after the inclusion of the emission line at £ ~ 6.4 
keV, the xl improves from 1.16 for 2242 degrees of freedom 
(dof)to 1.08, for 2239 dof(F-test of 9.2x10"^^). Although resid- 
uals improve significantly above 3 keV (see Fig.|4]third panel), 
other emission features can be seen at ~2.4, ~1.9 keV and ~1.3 
keV. Including an emission line at ~2.4 keV improves the fit fur- 
ther resulting in a xl of 1.07 for 2236 dof (F-test of 9.5 x 10"^ 
see Fig. |4] fourth panel). Adding another one at ~1.9 keV, we 
obtained a xl of 1.06 for 2233 dof (F-test of 2.6 x lO"*", see 
Fig. ID fifth panel). Finally, the last emission line at ~1.3 keV 
marginally improves the fit, leading to xl = 1 -05 for 2230 dof 
(F-test of 0.4, see Fig.|4]bottom panel). 

Using the combined spectra, we did not obtain an im- 
provement in the fit quaUty (xl - 1 .27, compared to xl - 
1.16 obtained using only the PN spectrum). Nevertheless, as 
evident from FigiH the emission lines are clearly seen in the 
raw data of the three cameras. In Table |3] we show the best 
fit parameters for the emission fines we detected in Fig. HI 
where the top spectrum in the top panel is from PN and the 
bottom spectra in the top panel are from MOS-I (black filled 
square) and MOS-2 (grey cross), respectively. 

The F-test is known to be problematic when used to 
test the significance of an additional spectral feature (see 
IProtassov & van Dik 2002l l. even if systematic uncertainties are 
not an issue. However, the low false alarm probabilities may 
make the detection of the line stable against even crude mistakes 
in the computation of the significance (Kxeykenbohm 2004|). 
Therefore, taking into account these caveats, we can conclude 




1 2 
channel energy (keV) 

Fig. 4. Fluorescence emission lines modeling of the phase- 
averaged spectra PN & MOS of 4U 1538-52. Top panel: Spectra 
and best fit model (three absorbed power-laws and four emission 
lines) obtained with EPIC/PN and EPIC/MOS. Bottom panels 
show the residuals in units of cr for models taking different num- 
bers of emission lines into account (MOS-1 black filled square 
and MOS-2 grey cross). Second panel: Without emission lines. 
Third panel: Fluorescence iron emission line at ~6.4 keV. Fourth 
panel: Fluorescence iron emission line plus a recombination line 
at ~2.42 keV. Fifth panel: Fluorescence iron emission line plus 
two recombination lines at ~2.42 keV and ~1.90 keV. Bottom 
panel: Fluorescence iron emission line plus three recombination 
Unes at -2.42 keV, -1.90 keV and -1.34 keV. 
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Table 3. Fitted parameters for the emission lines detected in 
Fig. fusing PN and MOS data simultaneously. 



Component 


Parameter 


3PL 


2PL+bb 


Fluorescence 


Energy (keV) 


6 4159+11.0011 

"•^^"'-0.008 


41 88+O-OOOt. 


iron line 


o-(eV) 


< 23 


<60 




C VV ^ C V ^ 


53 ±5 


52 ±5 


Identification 




Fe Ka I- 


-XVII 




F-test 


9.2x10-3^ 




Other 


Energy (keV) 


9 499+0.0023 
^•^^^-0.0020 


9 49 -J +0.0023 
^■^^■-'-0.0019 


emission lines 


o-(eV) 




C+50 




EW (eV) 


26 + 8 


27 ± 9 


Identification 




S XV Hea? 




F-test 


9.5x10-" 






Energy (keV) 


1 QfAj: +0.0024 


1 Qc 1 +0.0020 




cr (eV) 




<40 




FW (pV^ 

J_i VV V J 


50+ 12 


40+8 


Identification 




Si XIII 


Hea? 




F-test 


2.6x10^** 






Energy (keV) 


1 149+0.015 


1 340+"'"'' 




o-(eV) 


< 21 


< 21 




EW (eV) 






Identification 




Mg Key, Mg 


XI Hea? 




F-test 


0.4 








1.27(3725) 


1.27(3725) 



that the fluorescence emission lines are detected with high sig- 
nificance in the spectrum of this source and none of these fea- 
tures is the result of poor calibration. 

However, none of these models deals with the strong feature 
at ~2.1 keV. We added a Gaussian absorption profile at this en- 
ergy and found an improvement of the fit, but the high optical 
depth and its unconstrained value prevent us from accepting it. 
We checked the possibility that this feature could be due to a gain 
effect by changing the gain and offset values by up to +3% and 
+0.05 keV, respectively. The line profile expressed as ratio of 
a spectral data to a best-fit model was always around 70%, sug- 
gesting that it is not a gain effect and the absorption feature could 
be real. Moreover, although the feature is less pronounced 
in a spectrum extracted using single pixel events only, the 
EPIC/PN calibration team have confirmed (M. Guainazzi, 
private communication) the reality of the absorption feature 
at 2.1 keV. 

4.3. Pulse phase resolved spectra 

The hardness ratio shown in Fig. |2] demonstrates variation 
of the spectrum with respect to the spin phase of the neutron 
star. We divided the pulse into ten phase intervals of equal 
length in order to investigate the X-ray spectrum as func- 
tion of pulse phase. We produced energy spectra correspond- 
ing to each of these phase intervals, which were fitted by the 
same X-ray underlying model used to fit the pulse phase av- 
eraged spectrum, i.e., three absorption power laws plus an 
iron emission line. The resulting fit parameters are reported 
in Table II 



Table 4. Results of the fit of the pulse phase resolved spectra in the energy range 0.3-1 1 .5 keV. 



Pulse phase 



p a 1- a m p tf^T* 


0.05 


0.15 


0.25 


0.35 


0.45 


0.55 


0.65 


0.75 


0.85 


0.95 


Hard component 
N„ (10^2 cm-2) 

I. llVJWJll. LLl\^\^A. 

Normalization" 


49 ±9 

93+012 

57_ii 


1 114+0.04 
on +40 


"'-10 

1 CI +0.09 

i-'i-o.n 

1 An+60 
160_,n 


47!^ 
1 65 + 12 

9in+90 


67+" 
"'-11 

1.19 ± 0.10 

in-? +23 
103_3o 


79!;3 

95+0 05 

"•^-^-0.09 


59+" 
79+0 03 

"•'^-0.09 

^;4+22 
64_i7 


48+' 
1 l6+oo« 
1 9+T 


58 ±8 

1 44+0.08 

■^•^-0.11 

1 sn+90 


62!^ 
1.31 ± 0.12 


Scattered component 
iV// cm ) 
Normalization" 


1 T n+2-3 

iZ.U 2 

21. 6«' 


1 A+1-2 
i J.O , 2 

44-^f 


1 9 7+1.4 
12.7 , 3 

54+* 


Q 9+0.8 

47+fS 
^'-14 


1 1.'+ ± l.U 

61!jo 


19 8-1-1 1 
IZ.O ±1.1 

49+10 


1 1 ^+0.8 
11.6. 5 

32+^ 


1 n s+0.7 
10.8, 4 

48+13 


19 n+i-2 
59+' 


10.4_,2 
34^ 


Soft component 

Nu CI cm~^ 
li nil ) 

Normalization" 


2 0+°'* 
^•"-n.i8 


76+" " 

"•'"-0.08 
9 CI +0.16 
^•-'^-0.21 


78+°'' 
"•'°-p.ii 
9 0+0.5 


Q 7^+0.12 
9 Q+0.7 


67+° '^ 

9 Cn+0.4 


70+011 

"•'"-0.09 

9 c+0.4 

^•-"-0.3 


'59+010 

91+-Sf' 
^•l-n.3 


69+012 

9 /:+0.5 
^•0-0.4 


90 + 1 2 

T £+0.6 
^•^•-05 


67+° '° 

"•"'-0.09 
9 /I+0.5 
^•4-0.4 


Iron line 






















Energy (keV) 
o-(eV) 


<70 


6 40+°"'* 

-fl>03 

30+'* 


6.37 ± 0.03 
<60 


TO+0.09 

150 + 60 


6.44 ± 0.05 
60 ±60 


6.44 ±0.11 

130+'°° 


6 39+0 03 
< 70 


6.39 ± 0.03 
< 50 


6 41+°°' 
70+* 

'"-30 


6 37+010 

loo^f 


EW (eV) 
Unabsorbed flux* 


54 ± 16 
1.63 


55+2° 

-^-^-16 

2.08 


59+2' 
1.95 


70 ±30 
2.06 


49+18 
2.27 


00_3Q 

3.05 


35+" 
2.48 


^"-12 
2.30 


67+20 

"'-19 

2.37 


60 ±30 
1.87 


C-statistic(dof) 


2362(2239) 


2289(2239) 


2367(2239) 


2426(2239) 


2337(2239) 


2337(2239) 


2255(2239) 


2302(2239) 


2389(2239) 


2362(2239) 


"(x 10 in units of photons s ^ cm ^ 


keV-' at 1 keV 


















''(X 10-'°) erg S-' cm-2 
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We did not find statistically significant residuals in the 
energy range between 1 and 3 keV, and therefore no low en- 
ergy emission fines are required to fit these spectra, although 
this could be because of the lower statistics of the phase re- 
solved spectra with respect to the phase averaged spectrum. 
The energy of the iron emission line is compatible with being 
unchanged along the pulse profile, although the variations in 
the depth, width and equivalent width could be due to the 
presence of other unresolved emission lines. The BeppoSAX 
spectrum showed similar results with the pulse phase of 
the iron emission line parameters using four phase intervals 
(Robba et al. 2001 1. In other pulse phase spectroscopy analy- 
sis either the iron emission line was not present significantly 
(EXOSAT; Robba et al. 1992) or was also unchanged when 
the pulse phase profile was divided into two phase intervals 
{Tenma; Makishima et al. 1987). The variation of the differ- 
ent column absorption parameters are significant along the 
pulse profile, but again these do not show a clear relationship 
to the other parameters. 







1 1 1 ■ 1 ■ 






+4-' 


: ^ + + ^ 




+ + ^ 








• , , 








H — 1 — 1 — 1- 





I 1 I I I I 1 I 1 I L_ 

0.5 1 1.5 



Pulse phase 

Fig. 5. Photon index (top panel) and the corresponding unab- 
sorbed flux in the 0.3-11.5 energy band (10"'" erg s"' cmT^, 
lower panel) versus the corresponding phase interval. 



The photon index presents a clear modulation with the 
pulse phase showing two maximum values at phases 0.35 
and 0.85 and one minimum value at phase 0.65. In Fig. |5] 
we plot the photon index (top panel) versus the correspond- 
ing phase interval together with the corresponding unab- 
sorbed flux in the 0.3-11.5 keV energy band (lower panel). 
As can be seen from this figure, the phase of maximum un- 
absorbed flux is close to the phase of minimum photon in- 
dex, however there is not a one-to-one between these two pa- 
rameters. The unabsorbed flux presents two flat regions in 
the phase intervals 0.15-0.35 and 0.65-0.85 and increasing 
or decreasing flux in the other interval phases, it peaks at 
phase 0.55 and is near the minimum in the photon index. The 
result obtained from BeppoSAX showed a clear anticorrela- 
tion between the photon index and the flux, with the spec- 
trum being flatter in correspondence with higher flux levels 
(Robba et al. 2001). Moreover, this behaviour of the photon 
index fro m BeppoSAX was s imilar to p revious results from 
EXOSAT dRobba et al. 199211 and Ginga dClark et al. 1990» . 



Table 5. Fitted parameters for the X-ray continuum of the eclipse 
spectra in Fig. fusing PN and MOS data simultaneously. 



Component 


Parameter 


2PL 


Scattered 


Nh (10^- cm-^) 
Normalization^ 


79 ±3 
84+8 


Soft 


Nh (10^' cm-2) 
Normalization" 






Photon index^ 
C-statistic(dof) 


1 74+0.04 
6452(3727) 



" (x 10 **) in units of photons s ' cm ^ keV ' at 1 keV 
Both power-laws with the same photon index 



4.4. Detecting emission lines in tine eciipse spectrum 

As we showed in Sect. |3] the eclipse of the X-ray source is 
clearly visible in the last 55 ks of the XMM-Newton observa- 
tion (see Fig. [T]). Therefore, we analyzed the phase-averaged 
spectrum out of eclipse and in eclipse. The X-ray continuum 
spectrum out of the eclipse could be described by the models 
used in this work. While the parameters of the fluorescent 
iron emission fine at 6.4 keV were consistent with those of 
the whole observation, we did not find evidence of the other 
emission fines. Nevertheless, Fig.[T]suggests the presence of a 
variabiUty out of the eclipse since the HR changes with time. 
Therefore, we divided the data out of the ecUpse into two 
temporal intervals, 0-10 ks and 10-20 ks, taking into account 
the HR variability. We could describe these spectra with the 
same X-ray continuum, i.e. three absorbed power-law com- 
ponents, and a fluorescence iron emission fine. No statisti- 
cally significant residuals were observed in the energy range 
between 1 and 3 keV, and therefore no more emission lines 
are required to fit these spectra. The photon index and the 
iron emission line were compatible within the associated er- 
rors, but the photoelectric absorption were significantly dif- 
ferent. 

The continuum spectrum in ecUpse was described by two 
absorbed power-laws due to the fact that the direct compo- 
nent should not be seen in the ecUpse. In Fig. |6] we show 
the eclipse spectra (top panel) and the residuals to the two 
absorbed power-laws continuum model from PN and MOS 
cameras simultaneously (second panel). Although the fit 
quaUty is very good, the high value deduced for the absorp- 
tion column is at odds with the interpretation of this power 
law as the scattered component. In fact, the values are more 
consistent with those of the hard component. However, any 
other combination we tried, resulted in significantly worse 
fits. Whatever the interpretation might be, this component 
seems to be the dominant one both in and out of ecUpse. 

4U 1538-52 shows some emission lines in the eclipse, when 
the X-ray continuum emission is at a minimum. Evident are the 
presence of the Fe Ka line at 6.4 keV, an absorption edge around 

7.1 keV and emission lines between 1 and 3 keV. We fixed the 
two absorbed power-law components and added Gaussian 
profiles to fit these emission lines. Table|5]Usts the parameters 
of the continuum model we used to fit the spectra. 

Fig.|6]shows the spectra and the residuals with respect to 
X-ray continuum model in the whole energy range 0.3-11.5 
keV using the PN and MOS cameras. We note that in the 6.4- 

7.2 keV energy band some iron emission fines are present 
in the residuals. As has been shown using high resolu- 
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Table 6. Fitted parameters for the iron emission lines and for the 
recombination emission lines of He- and H-like species detected 
in the eclipse spectra (Fig.|6|l. 



Component 



Parameter 



Fluorescence 
iron line 

Identification 

Recombination 
emission lines 

Identification 



Identification 



Identification 



Identification 



Identification 



Energy (keV) 
o-(eV) 
EW (eV) 



Energy (keV) 
a-(eV) 
EW (eV) 



Energy (keV) 
o-(eV) 
EW (eV) 



Energy (keV) 
o-(eV) 
EW (eV) 



Energy (keV) 
o-(eV) 
EW (eV) 



Energy (keV) 
cr(eV) 
EW (eV) 



C-statistic(dof) 



6.4050 



2PL 

+0.005 
0.0010 

< 10 



262+21 
Fe Ka I-XVII 



6 634+° "-° 

"•"-'^-0.010 

<49 

84+12 

FeXXV 



? 4477+0.0023 

^•^^^'-0.023 

< 18 

47+15 

S XV Heor 



2 00 r''-''^^ 

■^•""^-0.05 

<280 

0+14 
°-8 

Si XIV Lyo- 



< 103 
21+" 
Si XIII Heff 



1.34 (fixed) 
30 (fixed) 
< 14 

Mg Ka, Mg XI Hett 
5148(3373) 



tion Chandra/HETG data (Brandt & Schulz 2000: D' Ai etaLl 
2007), the Fe complex can be clearly resolved into the near 
neutral Fe Ka line plus hot lines from highly ionized species 
of Fe XXV and Fe XXVI. Therefore, we used Gaussian lines 
to fit these residuals. In the bottom panels of Fig.|6]we show 
the ratio before and after adding the emission fines described 
in this Section. In Table |6] we report the parameters of the iron 
emission lines. The parameters of the Gaussian emission lines 
are Energy, cr and EW, indicating the centroid, the width and the 
equivalent width, respectively. 

The EW of the neutral iron Ka line is noticeably 
larger in eclipse (262 eV) than out of eclipse (~30 eV). 
|Rodes et al. (2006) found that the EW was larger when the 
source flux was low by using an RXTE observation which cov- 
ers nearly a complete orbital period. The measured line energy 
is consistent with an ionization stage up to Fe XVII, and this re- 
quires the ionization parameter ^ should be less than some hun- 
dreds ( K allman & McCray 1982t lEbisawa et al. 1996 ) . 

On the other hand, the simultaneous presence of the Fe XXV 
line in the spectrum implies an ionization parameter for the pho- 
toionized plasma of ^ ~ lO^*^ erg cm s"' (Ebis awa et al. 1996b . 
The broad range of the ionization parameter suggests either 
that the emitting material is present over a wide range of 
distances from the neutron star or has a large range of den- 
sities. 

Fig.|6]also shows the residuals with respect to a two absorbed 
power-law continuum model in the 1.0-3.0 keV energy range. 
Discrete recombination lines from He- and H-like species are 



present in the MOS and PN spectra of 4U 1538-52. The energy 
resolution of the instruments are not sufficient to resolve the lines 
clearly (Drake 1988; Hojnacki et al. 2001). Therefore, our iden- 
tification could be a blend of other emission lines. We also fit 
these emission lines as Gaussians and list them in Table |6l 

5. Summary and discussion 

We presented the spectral analysis of the HMXB 4U 1538-52 
using an XMM-Newton observation. The X-ray continuum is 
well fitted by three absorbed power-laws with a photon index 
~1.13, describing the hard, scattered and soft excess, respec- 
tively. The infeiTed unabsorbed flux is ~ 2. 1 x 10"'" erg s"' cm"^ 
in the 0.3-1 1.5 keV energy band, coiTesponding to a luminosity 
of ~ 7.5 X 10^^ erg s"', assuming an isotropic emission and a 
distance to the source of 5.5 kpc (B ecker etal. 1977) . The flux 
found by RXTE in the 3-11.5 keV energy band and an orbital 



phase of 0.85 was 1.3x10" erg s" cm" . Using the spectrum 
obtained by XMM-Newton over the same range of phases, the 
flux we obtained in this work was two times lower, 6.3x10"" 
erg s"' cm"^. 

The soft excess is present in the spectrum and can be mod- 
elled with different absorbed components. We simply modelled 
the soft emission with an absorbed power law component, al- 
though Fig. ID still showed spectral residuals at lowest energies. 
Our results show that a blackbody component could also 
be the physical origin of this soft excess, taking the asso- 
ciated errors into account. The soft excess in other HMXBs 
has been explained by a blend of Gaussian emission lines only 
(^Boroson et al. 2003b . We also tried to fit the soft emission using 
Gaussian profiles, but we did not obtain a significant improve- 
ment of the fit because of the low level of counts below 0.6 keV. 

We detected an iron Ka line at ~6.41 keV, with an EW of 
~50 eV. The BeppoSAX observation of this system obtained 
an EW of 57 eV in the same orbital phase range 0.75-1.00 
dRobba et al. 200 il l and showed an increase in the post-egress 
phase to 85 eV. The phase-averaged spectrum obtained by RXTE 
reported an EW of 62 eV (Coburn et al. 2002). In addition this 
iron line is detected in all orbital phases (Robba et al. 200T1 
Rod es 2007b . therefore the Ka iron line is not only produced by 
fluorescence from less ionized iron near the neutron star's sur- 
face but also a fraction of the observed line flux must originate 
from more extended regions. We have also detected a num- 
ber of emission lines which we interpret as recombination 
fines from highly ionized species. Since these fines are de- 
tected in ecfipse, they must be produced in an extended halo. 
Likewise, we have found an absorption feature at 2.1 keV. 
Whether it is produced by physical properties of the source 
or it is due to cafibration elfects, is stUl an open issue. 

We compared the phase-average spectrum to the eclipse 
spectrum. In the phase-average spectrum, we found no evidence 
of any other iron line apart from that at 6.4 keV, and the absorp- 
tion edge at ~7.1 keV was well described by the X-ray absorp- 
tion model. The 6.4-7.2 keV energy band showed a complex 
structure, but we did not find a proper model to describe it or 
detect other iron features significantly. We also detected discrete 
recombination lines in both EPIC/PN and EPIC/MOS spectra. 
The emission lines reported in Sect. 14.21 can be identified with 
He- and H-like ions. Although we could identify some of them, 
the energy resolution of the EPIC instruments is such that the 
lines could be blended. 

However, in the eclipse spectrum, we detected the fluores- 
cence iron line at 6.4 keV and one more iron emission line at 
6.6 keV. The discrete recombination lines detected are listed in 
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Fig. 6. Iron emission lines and recombination emission lines of He- and H-like species in the eclipse spectrum of 4U 1538-52. 
Top panel: Spectra and best fit model (two absorbed power-laws and six emission lines) in the 0. 3-11. 5 keV energy range obtained 
with PN (top) and MOS (light grey filled squared and dark grey cross) cameras. Bottom panels show the residuals. Second panel: 
Without emission lines. Third panel: With iron and recombination of He- and H-like emission lines for the combined spectra. 



Table |6] The presence of these lines in an eclipse spectrum im- 
plies that the formation region extends beyond the size of the 
B supergiant. Moreover, the ionization state was estimated to 
range from 10^ ' to lO-''^ erg cm s"', due to the simultaneous 
detection of elements with both low and high ionization levels. 



This broad range of ^ also suggests either that the emitting 
material is present over a wide range of distances from the 
compact object or has a large range of densities. 

The pulse phase-resolved spectroscopy showed signifi- 
cant variability of the photon index and the unabsorbed flux, 
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but no clear correlation or anticorrelation between them. 
Significant variations with the pulse phase were also ob- 
served in the difTerent column absorption values, but again 
did not show a clear relationship to the other parameters. 

Future observations with high spectral resolution instru- 
ments will be needed to unambiguously resolve possible 
blended lines found in this study allowing the full use of their 
diagnostic capabiUty. 

Acknowledgements. Part of this work was supported by the Spanish Ministry of 
Education and Science Primera ciencia con el GTC: La astronomia espanola en 
vanguardia de la astronomia europea CSD200670 and Multiplicidad y evolucion 
de estrellas masivas project number AYA200806166C()303. This research has 
made use of data obtained through the XMM-Newton Science Archive (XSA). 
provided by European Space Agency (ESA). We would like to thank the XMM 
helpdesk, particularly Matteo Guainazzi, for invaluable assistance in determin- 
ing the systematic uncertainties in the PN data. KLP and JPO acknowledge sup- 
port from STFC. JMT acknowledges the support by the Spanish Ministerio de 
Educacion y Ciencia (MEC) under grant PR2007-0176. JJRR acknowledges the 
support by the Spanish MEC under grant PR2009-0455. 



References 

Arnaud, K. A. 1996, in Astronomical Data Analysis Software and Systems V, 

ed. J. H. Jacoby & J. Barnes, ASP Conf. Ser. 101, San Francisco, 17 
Audley, M. D., Kelley, R. L., Boldt, E. A. et al. 1996, ApJ, 457, 397 
Becker, R. H., Swank, J. H., Boldt, E. A. et al. 1977, ApJ, 216, LI 1 
Brandt. W. N. & Schulz, N. S. 2000, ApJ, 544, LI 23 
Boroson, B., Vrtilek, S. D., Kallman, T. et al. 2003, ApJ, 592, 516 
Clark, G. W., Woo, J. W., Nagase, F. et al. 1990, ApJ, 353, 274 
Clark, G. W., Woo, J. W. & Nagase, F 1994, ApJ, 422, 336 
Clark, G. W. 2000, ApJ, 542, L131 

Coburn, W., Heindl, W. A., Rothschild, R. E., et al. 2002, ApJ, 580, 394 

Corbet, R. H. D., Woo, J. W. & Nagase, F 1993, A&A, 276, 52 

D'Ai, A., laria, R.. Di Salvo, T. et al. 2007, ApJ, 671, 2006 

Davison, R J. N., Watson, M. G. & Pye, J. P 1977, MNRAS, 181, 73P 

Drake, G. W. 1988, Canadian J. Phys., 66, 586 

Ebisawa, K., Day, C. S. R., Kallman. T. R. et al. 1996, PASJ, 48, 425 

Eisner, R. F & Lamb, F K. 1977, ApJ, 215, 897 

Giacconi, R., Gursky, H., Kellog, E. et al. 1974, ApJ Suppl., 27, 37 

Harding, A. K. & Daugherty, J. K. 1991, ApJ, 374, 687. 

Hatchett, S. & McCray, R. 1977, ApJ, 211, 552 

Hickox, R. C, Nai-ayan, R. & Kallman, T. R. 2004, ApJ, 614, 881 

Hojnacki, S. M., Kastner, J. H., Micela, G. et al. 2007, ApJ, 659, 585 

laria, R., Di Salvo, T., Robba, N. R. et al 2005, ApJ, 634, L161 

Kallman, T. R. & McCray, R. 1982, ApJS, 50, 263 

Kreykenbohm, 1. 2004, Ph. D. thesis. University of Tubingen 

Makishima, K., Koyama, K., Hayakawa, S. et al. 1987, ApJ, 314, 619 

Mon-ison, R. & McCammon, D. 1983, ApJ, 270, 119 

Mukherjee, U., Raichur, H., Paul, B. et al. 2006, JAA, 27, 41 1 

Nagase, F, Corbet, R. H. D., Day, C. S. R. et al. 1992, ApJ, 396, 147 

Protassov, R. & van Dik, D. A. 2002, ApJ, 571, 545 

Reynolds, A. R. Bell, S. A. & Hilditch, R. W. 1992, MNRAS, 256, 631 

Robba, N. R., Cusumano, G., Orlandini, M. et al. 1992, ApJ, 401, 685 

Robba, N. R., Burderi, L., Di Salvo, T. et al. 2001, ApJ, 526, 950 

Rodes, J. J. 2007, Ph.D. thesis. University of Alicante, 

http://hdl.handle.net/10045/13227 
Rodes, J. J., ToiTejon, J. M. & Bernabeu, G. 2006, Proceedings of the The X-ray 
Universe 2005, 26-30 September 2005, El Escorial, Madrid, Spain. Editor: 
A. Wilson (Noordwijk: ESA Publications Division), ESA SP-604, 1, 287 
Rodes, J. J., Torrejon, J. M. & Bernabeu, G. 2006, The X- 
ray Universe 2008, 27-30 May 2008, Granada, Spain, 
rhttp://xmm.esac.esa.int/extemal/xmm_science/workshops/20Q8symposium/# topicB | 
Rodes-Roca, J. J., Toirejon, J. M., Kreykenbohm, I. et al. 2009, A&A, 508, 395 
Sako M., Liedahl D. A., Kahn S. M. et al. 1999, ApJ, 525, 921 
Struder, L., Briel U., Dennerl, K. et al. 2001, A&A, 365, L18 
Turner, M. J. L., Abbey, A., Ar-naud, M. et al. 2001, A&A, 365, L27 
van der Meer, A., Kaper, L., Di Salvo, T. et al. 2005, A&A, 432, 999 
van Loon, J. Th., Kaper, L. & Hammerschlag-Hensberge, G. 2001, A&A, 375, 
498 

White, N. E., Swank, J. H. & Holt, S. S. 1983, ApJ, 270, 711 



